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The overlap fermion offers the tremendous advantage of exact chiral symmetry on the lattice, but 
is numerically intensive. This can be made affordable while still providing large lattice volumes, 
by using coarse lattice spacing, given that good scaling and localization properties are estabhshed. 
Here, using overlap fermions on quenched Iwasaki gauge configurations, we demonstrate directly 
that the overlap Dirac operator's range is comfortably small in lattice units for each of the lattice 
spacings 0.20 fm, 0.17 fm, and 0.13 fm (and scales to zero in physical units in the continuum 
Umit). In particular, our direct results contradict recent speculation that an inverse lattice spacing 
of 1 GeV is too low to have satisfactory localization. Furthermore, hadronic masses (available on 
the two coarser lattices) scale very well. 
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1. Locality 

In the last several years the use of the overlap fermion has become more popular because the 
conceptual and technical clarity that results from its exact chiral symmetry on the lattice is seen 
to overcome its superficially higher computational cost as compared to conventional formulations 
of the fermion. Furthermore, the disparity in numerical intensity can be mitigated by using coarse 
lattice spacing, given that good scaling and localization properties ai^e established. 

Hernandez, Jansen, and Liischer ^ showed numerically that Neuberger's overlap operator is 
local (using the Wilson gauge action on fine lattices). Recently, Golterman, Shamir, and Svetit- 
sky 1^ have speculated that overlap simulations with a cutoff of IGeV (such as [^]) might have 
a range as long as 4 lattice units, and thus be afflicted by unphysical degrees of freedom as light 
as 0.25 GeV. Here we show directly that such is not the case; the range is about 1 lattice unit (in 
Euclidean distance or 2 units of "taxi-driver" distance). All is well. 

1.1 Lattice Details 

We use the renormalization-group-improved Iwasaki |Q] gauge action, on three different lat- 
tices; for each, the lattice size, lattice spacing, and number of configurations used are tabulated in 
Table |l|. 

A';, X Nt a(fm) N^fg 



16^x28 0.20 300/10 
20^x 32 0.17 98/10 
28^x44 0.13 /lO 



Table 1: Lattice size, lattice spacing, number of configurations (for scaUng/locaUty). 



For the associated scaling study of hadron masses, we use the overlap fermion ^ and massive 
overlap operator |^ 

^K) = (p + ^) + (p-^)75s(//) 

where £{H) = H /VlP-, H = YsD^, and D„ is the usual Wilson fermion operator, except with 
a negative mass parameter — p = 1/2k" — 4 in which < K < 0.25; we take K = 0.19 in our 
calculation which corresponds to p = 1.368. For the locality study, we set mo = to look at the 
properties of the massless operator D{0). Complete details are described in [^. 

1.2 Locality as Measured by Taxi-Driver Distance 

It is convenient for formal reasons to discuss locality in terms of "taxi-driver" distance 

rTD = \\x-y\\i =Y,\x^-y^\ (1.1) 
J" 

The locality of the overlap operator is then studied by plotting the quantity \D{r)\ (f{r) in the 
notation of |[l]]) as a function of the taxi-driver distance for a localized source, Ya{x) = 5{x)5ap 
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for fixed Dirac-color index j8 . 

|D(r)| =max{||Di/^(x)|| | = r} (1.2) 

For large r, the kernel of the Dirac-overlap operator decays exponentially with decay rate v = r^^ , 
where j-qv is the range (characteristic decay distance) measured in lattice units. 

1.3 Results 

In the left pane of Fig. |l|, we plot \D{r)\ as a function of taxi-driver distance for each of three 
lattice spacings. At large distances, we fit to an exponentially decreasing function to extract the 
range Toy. These are tabulated in Table ^. We note that our results are consistent with the results of 
Hernandez, Jansen and Liischer [|I]] on finer lattices for the overlap operator (p = 1 .4) with Wilson 
action at j8 = 6.0, 6.2, and 6.4 where they find V = r^^^ = 0.49. 

0.20fm 1.93(1) 
O.lVfm 1.83(1) 
0.13fm 1.81(1) 

Table 2: The range (taxi driver metric) for three lattice spacings. It is less than two lattice units on a lattice 
as coarse as 0.20fm. 
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Figure 1: Left: For each of three lattice spacings, the expectation value of \D{r)\ as a function of the 
taxi-driver distance, r. For large r, \D{r))\ falls exponentially, with range Toy. Fitted values of rov are shown 
with fit intervals. Right: Taxi driver range in physical units (fm) as a function of lattice spacing. The range 
is small even at coarse lattice spacing and trends to zero in the continuum limit. 

Furthermore, it is gratifying to see that even for our coai^sest lattice, 0.20fm, the measured 
range is less than two lattice units. In the right pane of Fig. |l] we plot the range in physical units as 
a function of lattice spacing; it trends to zero in the continuum limit. 
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1.4 Locality as Measured by Euclidean Distance 

We conclude that it is perfectly acceptable to simulate overlap fermions with lattice spacing 
as coarse as 0.20 fm, since for this we find that the range is not greater than two lattice units when 
measured in taxi-driver distance. In fact, the situation is even better than it seems. To see this, we 
consider the more familiar standard Euclidean metric 



ni = \\x-y\\2= n^\x^-yii\'^ (1.3) 
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Figure 2: For the 16^ x 28 lattice, the expectation value of \D{r) \ as a function of Euclidean distance r (for 
data cut to remove wrap-around effects). Data is plotted at each Euclidean distance in red. The blue points 
are averages within bins of width Ar = 0.5. The exponential tail is then fit over the longest interval possible 
with acceptable shown by the green Une which has inverse slope rov = 1 .05 ( 1 ), the range. 

In Fig. ^ we plot D{r) versus the Euclidean distance for our coarsest lattice. As compared to 
the corresponding plot [l] using the less-famihar taxi-driver distance, one sees that the data are more 
scattered due to violations of rotational symmetry, but are still clearly contained with a worst-case 
decay rate. 

0.20fm 1.05(1) 
O.lVfm 0.98(1) 
0.13fm 0.9(1) 

Table 3: The range (Euclidean metric) at three values of lattice spacing. It is less than about 1 lattice unit 
with lattice spacing as coarse as 0.20 fm. 
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Again we fit tlie tail of the function with a decaying exponential to extract the range. We 
tabulate the results in Table ^. Note that although the two ranges (taxi-driver and Euclidean) differ 
by a factor of two they are quite compatible heuristically; on a L"* hypercube, the maximum taxi- 
driver distance is 4L, and the maximum Euclidean distance is \/4Z? = IL. 

So even at lattices as coarse as a = 0.2 fm, the range is about 1 lattice unit (measured using 
Euclidean distance, or 2 units using taxi-driver distance). No unphysical degrees of freedom are 
induced at distances longer than the lattice cutoff. 

2. Scaling 

At Lattice 2004, Davies et al. collected world data to demonstrate that different quenched 
quark formulations could have a consistent continuum limit. 

quenched comparison 
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Figure 3: "Aoki" plot for various quenched data, as obtained from [^. Our data is labeled "overlap". 

The conclusion, as illustrated in Fig. |3|, is that they could. But we emphasize that the con- 
strained fit demanded that there exist a global continuum limit (by design). Furthermore, the global 
continuum limit differed substantially from the continuum limit obtained solely from the Wilson 
formulation where the discretization eiTors are largest. The lesson learned is that with large dis- 
cretization errors, it is quite possible to be misled when extrapolating to the continuum limit, even 
with high statistics and many lattice spacings. It is important to seek a formulation with very small 
discretization errors to be able to trust the continuum extrapolation. 
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Here we add our data "overlap" to the global quenched spectrum data and find that of all the 
formulations, its discretization eiTors are smallest, allowing for viable computation at surprisingly 
coarse lattice spacing. 

As another example of the efficacy of the overlap formulation, we made a non-perturbative 
computation ^ of the renormalization constants of composite operators on the 16^ x 28 lattice 
using the regularization independent scheme. 

We found that the relations Za = Zy and Z5 = Zp agree well (within 1%) above m = 1.6GeV. 

The 

'WAqcd^'^ and (ma)'^ corrections of the renormalization are small; the mass dependence is less 
than about 3% up to ma = 0.6. This is much superior to the competition. 

3. Conclusions 

It is viable to simulate quenched overlap fermions at surprisingly coarse lattice spacing. Lo- 
cality is well under control; the range (characteristic exponential decay length) is about one lattice 
unit (of Euclidean distance, or about two lattice units of Taxi-driver distance) for lattice spacing as 
coarse as 0.20fm (such as in and trends to zero (in physical units) in the continuum limit. Scal- 
ing is remai^kable. The Aoki plot is essentially flat up to 0.20fm. The overlap fermion outperforms 
all other formulations; discretization errors are smallest for overlap. Non-perturbative renormaliza- 
tion of operators show little mass dependence [^; e.g. less than about 3% up to ma = 0.6 for the 
renormalization constants 
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